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INTRODUCTION 



Deardorff's (1972) similarity theory for convective boundary 
layers (CBL)s suggests that CBL parameters are controlled by 
buoyancy, B = g/0, (where g is gravitational acceleration and 0 
is potential temperature), surface temperature flux, w'0' o , 
height z above the surface, and z t , the height of the inversion 
base which marks the top of the mixing-layer (ML) . The velocity, 
temperature, and humidity scales: w,, u, ML , ©, N , q» ML , formed from 
these parameters are used in universal relationships for the 
gradients and fluxes of mean and turbulent quantities in the CBL. 
The first of these, the convective turbulence scaling velocity 



determines much of CBL physics, including diffusion. Along with 
w'G'o, w, forms the basis for the remaining parameters used to 
scale the heat, moisture, and momentum, as well as passive scalar 
fluxes. However, w, neglects mechanical turbulence due to shear 
and is only strictly valid for free convection induced by surface 
heating. Yet, w* is used often in mixed forced/free convection 
cases for lack of alternatives. Here we seek w, s , an amendment 
to w», which includes mechanical turbulence contributions under 
mixed forced/free convection conditions. We assume mechanical 
generation comes from three sources: the barotropic, horizontally 
homogeneous surface layer (SL), barotropic boundary layer 
entrainment zone (EZ) , and a baroclinity (BC) component. Beyond 
the z and z, length scales, we propose that mixed forced/free 
convection also depends on the SL roughness length, z 0 , SL 
Obukhov length, L = -u, 3 ©/ (gkw ' © ' 0 ) (where u, is friction velocity 
and k = 0.4 is the von Karmen constant), EZ thickness, Ah, and 
z b , a baroclinic length scale which we define below. We model 
these mechanical sources, using relevant ratios of the length 
scales, thereby fulfilling similarity requirements. 



We outline the basic approach before discussing details. For 
non-baroclinic , horizontally homogeneous, steady state cases, the 
vertical profiles of the buoyancy and momentum fluxes must be 
linear, if the shapes of the vertical profiles of potential 
temperature and wind speed are to remain static. That is, 



w. = (Bz, w'©' 0 ) 1/3 



(1) 



APPROACH 



Bw ' © ' z = Bw ' © ' 0 ( 1 - (1 + a,)z/Zj) 



(2a) 



-u,'w' z = -u,'w , 0 (l - (1 + a 2 ) z/z,) 



(2b) 



2 



where a, = -w ' 9 1 zi /w ' 8 ' 0 and a 2 = Uj ' w ' zi /Uj ' w ' 0 (see fig. 1). The 
subscripts "zi" and "0" refer to conditions at the inversion base 
and surface. -u'w' 0 = u, 2 , w'8' zi - WeA0, and -Uj'w' z ; = AU,We, 
where the CBL entrainment rate, We = z it and a9 is the potential 
temperature jump across the EZ. The subscript ",t" denotes the 
partial derivative with respect to time. Such notation will 
refer to spatial and temporal derivatives in general. The Betts 
(1974) correction, w'8' zi = (a0 - r u Ah)We, hardly affects eqn. 2a 
since Kamada (1988a) showed that the much smaller interfacial 
thickness, l m , must replace the EZ depth, Ah. For lack of better 
information, we assume that above z t the barotropic flux profiles 
are linear with height up to H = Zj + Ah, the top of the boundary 
layer . I . e . , 



z - Zj 

Bw'©' z = Bw ' 0 ' zi ( 1 - ) , (3a) 

Ah 

Z — 2 ■ 

-u,'w' z = -u.'w'jjU - ) . (3b) 

Ah 

Using eqns. (1), (2a), and (3a), the vertically integrated 

buoyancy flux or total CBL buoyancy production rate of turbulence 
kinetic energy (TKE) can be written to show that the purely free 
convective turbulent scaling velocity is proportional to the cube 
root of the boundary layer integrated TKE production rate, 

H 1 - a, ( 1 + a) 

j Bw'8' z dz = n w, 3 , where (1 = , 

z0 2 

(4) 



and a = Ah/Zj. n = 0 . 4 for mature CBLs because under such 
conditions, a, - 0.2 and a ~ 0.1. Since this may not be true for 
other conditions, we discuss n in detail below. However, here we 
propose that this sort of proportionality extends to shear; that 
mechanical contributions to a shear inclusive turbulent scaling 
velocity, w, s , may be written in terms of the boundary layer 
integrated mechanical generation rate of TKE. I.e., 



w 



*s 



$ Bw'8' z dz + j -Uj ' w ' u, Z dz 
zO zO 



(5) 



n 
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fig.l Schematic profiles of a) buoyancy, b) stress, c) U velocity, 
and d) TKE shear generation rate for barotropic case 



z z z 




We also propose that the CBL shear generation rate can be 
represented by discrete barotropic surface layer (SL) , 
entrainment zone (EZ) , and baroclinic (BC) components. 



EZ BC 

H H 

y -u, ’ w ' Ui z a Z + J-Uj'w'Ujb ^Z . 
zi 0 

( 6 ) 

Here u ib>z is just the baroclinic portion of the vertical shear. 



BAROTROPIC SURFACE LAYER SHEAR 

We first consider the barotropic SL and assume that most of the 
shear occurs below the height, -L. From SL similarity we have 

u z = (u,/kz)0 m , with <p m = (1 - 15z/L)' 13 , (7a, b) 



(Carl et al. 1973). Some workers prefer 0 m formulas involving 
the -1/4 power. However, the choice is among empirical 
expressions and integrals of -1/4 power expressions result in 
awkward complex numbers. The accepted value in the free 
convection limit is -1/3, so we choose this for theoretical 
appeal as well as mathematical convenience. 

Since <p m is only valid for z >> z 0 , we set the upper integration 
limit at -L, the lower limit at a somewhat arbitrary 10 z 0 , and 
treat the layer below 10 z 0 separately. Above 10 z 0 , from eqns. 
(2b) and (7a) , we have 



SL 



J— u,'w'u lz dz = f-Ui'w'u^dz + 
zO zO 



-L 

f -U/w'Uj flz 

lOzO 



U, 3 -L <p m dz 
\ 

k lOzO z 



u»(u , 2 - AUjWe) 



kZj 



J 



lOzO 
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U« 3 -L <P m dz -L 

\ - (l-a 2 ) (L/z i2 )w. 3 \ <p m dz 

k 10z0 Z 10z0 . (8) 



From the L definition, u, 3 /k = (-L/Z;)w» 3 . With the axes aligned 
with the mean wind, we have -u'jW'u., — -u'w'u z , since (v z , w z ) << 
u z in the mixed layer. If shear supplies sufficient TKE, 10 z 0 
<<< -L. So in eqn. (8) term (2) — 0.54 (1 - a 2 ) (L/Zj)^, 3 . 

Letting p = -15z/L, the first integral reduces to 



4 



- V3 ATAN 



V3 (2 (p+1) 1/3 + 1) 



Ln [ (p+l) 2/3 + (p+1) 1/3 + 1] 



Ln [ (p+l) ,/3 



1 ] 



-L 

lOzO 



The upper limit evaluates to = 1.51. The first two terms in the 
lower limit — 1.26, if L >> 150z 0 . So in eqn. (8) the first 
integral = 0.25 - Ln R, where R = (1 - 150z 0 /L) 1/3 - 1 << 1. A 
5cm z 0 and 150m L would give 150 z 0 = -0.05L. Above 10 z 0 , this 
leads to = 4.36u, 3 /k from term (1) in eqn. (8). Oddly, R varies 
inversely with roughness length. But from u, = ku/(ln(z/z 0 - y> m ) 
and the definition of L, we see that 



nz 0 /L = (ngw ' 0 ' /k 2 u 3 0 ) z 0 ( ln ( z / z 0 ) - <^ m ) 3 



(9) 



So for given levels of u, w'9', 0, and z, the ratio nz 0 /L and 
thus term (1) may vary much less than z 0 itself. 

Below 10 z 0 we neglect non-neutral contributions to u, and assume 
a constant value, -Uj'w'u z = u, 3 /10kz 0 . Then 

10z0 

f -u,'w'u z dz ~ u* 3 /k . (10) 

0 

Together, eqns. (8) and (10) give an estimate of the vertically 
integrated SL shear production rate, 



-L 

j -u,'w'u z 3z 
0 



(1.25 - Ln R) u» 3 /k , 



= (-w, 3 L/Zj) (1.25 - Ln R) . (lla,b) 

Note that we have dropped term (2) from eqn. (8), since 1) the 
coefficient is an order of magnitude smaller than for term (1), 
2) term (2) is second order in z ; , 3) usually, -L < 0.4 Z;, and 

4) a, is of order 10' 1 - 1; so some further cancellation occurs. 

From eqns. (2) and (3), if n = 0.4, we have the following ratio 
of enhanced to unenhanced TKE contributions, 



5 



n 



(1.25 - Ln R)L/Zj 



1/3 



n 



w, 



(1 + (3.1 - 2 . 5Ln R) (-L/Zj) ) l/3 w, . 

(12) 



For our previous example, w, s exceeds w, by a factor of — (1 + 

3 3u« 3 /w« 3 ) 1/3 . In this case, since u,/w, = 0.54, SL shear 
generation alone may enhance w, s by more than 80%. 

To be useful, w« s should maintain the same proportion to the mean 
turbulent velocity, a ui / 3, in both the convective and neutral 
limits. We test this by noting for the barotropic case that the 
cube root of the vertically integrated production rate is usually 
— 2.4 u, for SL shear. The EZ and BC terms are not relevant for 
the pure neutral case. From Grant (1986) and Mason and Thompson 
(1987), we find that a ui /3 = 1.3 u, in the layer averaged neutral 
boundary layer. This implies that a ui /3 — 0.53 w. s . In the pure 
convective limit, the CBL averaged value for a ui / 3 is " 0.55 w, 
(Caughey and Palmer, 1979 and Deardorff et al. 1980). So, the 
ratio of w» s to the layer averaged mean turbulent velocity is the 
same in the convective and neutral limits, as required. 



ENTRAINMENT ZONE SHEAR 

The EZ also generates mechanical turbulence. To quantify it we 
assume, moreover, that shear is constant across ah. Thus, 

H 

J -Ui'w'u iz dz = Weauf/2 . (13) 

zi 

Equation (13) suggests the need for a simple shear inclusive 
model to estimate the entrainment rate, We. For this purpose we 
parameterize the TKE equation, 



storage buoyancy transport 



shear dissipation 



e t = 6 w ' 9 ' + (u/e + p'e/p) z - Uj'Uj'Ujj -^ u i,j u i,j 



at the inversion base as, 



(14) 



storage buoyancy transport shear dissipation 

C st w, 2 We/Ah = -BA0We + C c w, 3 /Zj + C m AUfWe/Ah - C d w, 3 /z; 

(15) 
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The scaling arguments are that 1) again w'0' zi = -WeA0 and 2) 
large mixed layer eddies drive the TKE transport rate into the 
EZ. Thus TKE transport should vary with the TKE (oc w* 2 ) times 
the large eddy turnover time scale, Zj/w, . 3) The shear 
production rate is modeled as the momentum flux, AU,We, times the 
local EZ shear, AUj/Ah. 4) To maintain a steady state inertial 
subrange, the dissipation plus buoyancy destruction rates must 
balance the large eddy TKE injection rate due to transport and 
shear production. So in the neutral limit, dissipation should 
scale like transport with the change rate of large eddy TKE, 
w, 3 /Zj. Since dissipation occurs at small scales, unaffected by 
buoyancy anisotropy, we retain the neutral form. For simplicity 
we neglect the small TKE storage rate term. C m estimates vary 
somewhat, (Pollard et al. 1981; Kato and Phillips, 1967, Tennekes 
and Driedonks, 1981), but we see little reason for departures far 
from unity, so we assume C m ~ 1 . 0 . We discuss C c and C d below. 

Tennekes and Driedonks (1981) presented a similar model, except 
that here we scale the shear production locally with Ah instead 
of Zj. Solving for the entrainment rate gives, 



(C c - C d ) 

We = w./Zj . (16) 

J3a 0 - Auy/Ah 



Therefore, from eqn. (13) we have 



h (C c - C d ) Ah 

j -U/w'Uj Z 3z = w, 3 , (17) 

zi ' 2 ( R b - 1) Z, 



where the bulk Richardson number, 

(3a 0a h 

R b = — . (18) 

AU 2 

The EZ is transitionally turbulent by nature. Thus, we expect R h 
to fluctuate about some critical threshold turbulence value, R cb . 
R cb is much larger than the 1/4 value accepted for the stable SL 
because EZ TKE is mainly maintained by large eddy transport from 
below. So shear generation is only part of the TKE injection 
rate, and the dominant eddies are also much larger than in the 
transitionally turbulent SL. Thus, the EZ's large negative 
buoyancy will destroy much of the TKE during the eddy cascade 
before it can be dissipated at small scales. For these reasons 
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R cb may exceed unity in the EZ. 

Indeed, our simulations to steady state using such an entrainment 
model suggest a typical value of R cb = 1.3. From the results of 
Turner (1968), Deardorff et al. (1980), and many others, we 
suggest that C c - C d = 0.25 in the convective limit. So in eqn. 
(17) the EZ shear contribution = 0.5 aw, 3 , where a = Ah/Zj. 

To assume that a = 0.1 - 0.2 may often suffice. Then, a, = 0.2, 
fl = 0.4, and the barotropic EZ contribution to w, s is only a few 
percent. But a can exceed unity during rapid entrainment into a 
near-neutral layer remaining from the previous day's CBL (Nelson 
et al. 1989). During such times, entrainment models suggest that 
a, may exceed 40% (Kamada, 1988a). Then fl = 0.1, and EZ 
enhancement of w„ s may exceed 0.4. With SL shear this may cause 
w, s to double. In this case a better diagnostic or prognostic 
estimate for a is needed. For a diagnostic estimate, we begin by 
assuming R b = R cb and solve eqn. (17) for Ah. This leaves unknown 
aUj and a 6 , the inversion temperature jump. a© is hard to 
specify more precisely than - +/- 30% because z ; and H are not 
that highly resolved. Yet, we know that a© grows with the 
boundary layer growth rate, H t , into the layer aloft (with lapse 
rate r u ) , and shrinks with the mixed layer warming rate, 0 mt . 

So, using the simplest steady state assumption, 



a© t = 



r u H ,t - 



0 



m,t 



= 0 



(19) 



The mixed layer warming rate is given by the surface heating 
rate, plus the entrainment rate of warm air from aloft, i.e., 



0 



m,t 



(w'©' 0 - w'©' zi ) /z; 



( 20 ) 



So, if we approximate H t as We ( 1 + a) , we have 



a© = w'Q'o/We - r u Zi(l + a) 



( 21 ) 



Using eqns . (16) and (18) and combining constants, this results 

in two solutions for a. We show only the physical root, 



a 



2 . 3 aUj" 

F u G Zi 2 




+ 



4 



J 



1 



2 



( 22 ) 
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Though not explicit, eqn . (22) implicitly involves a© because it 

presumes through eqns. (18 - 22) that a0 depends on r u . We, 
w'0' o , and z t . This assumes that the capping inversion requires 
finite potential temperature and windspeed jumps, whenever there 
is horizontal wind. 

In eqn. (22) we must still estimate the windspeed jumps across 
the EZ. Wyngaard (1988) proposed the following first order 
estimate. For steady flow over a horizontally homogeneous 
surface, the stress gradient is given by 



-u'w' z = f ( v - v g ) , and -v'w' z = f(u g - u) , (23a, b) 

where f refers to Coriolis forcing, and the subscript, g, refers 
to geostropic winds (Panofsky and Dutton, 1984) . For barotropic 
flow with the axes aligned with the mean mixed-layer wind, 



-u'w' z = fv g = (We(u g - u) + u» 2 )/Zj , and 

-v'w' z = f(u g - u) = We ( v g - v ) / z j . (24a, b) 

The momentum balance becomes 

We(Ug - u) + u, 2 = -fz.Vg , and 

-We Vg = f z j ( u g - u) . (25a, b) 

For barotropic cases, if au = (u g - u) and av = (v g - v) across 
the EZ, then the solutions to eqns. (25a, b) are to first order, 

au = We(u,/fZj) 2 , and 



AV = U, 2 /fZj 



(26a, b) 



Here, almost all of the EZ shear results from turning the wind 
into the v direction. So for a mature mid-latitude BL with u, = 
0.4ms 1 , f = lO'V, z, = 10 3 m, 0 = 300 °C, and T u = 3 x 10' 3 °C m' 1 , 
we have a = 0.057. However, during growth into the near neutral 
remnant mentioned earlier, conditions may be more like r u = 1 x 
10" 3 °C m' 1 , and z K = 500m. This would result in a = 1.3. Values 
spanning such a range were reported in Nelson et al. (1989). 

They also observed and prognostically modeled the typical early 
morning to late afternoon hysteresis found in plots of a versus 
We/w,. They noted that such hysteresis is absent from previous 
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diagnostic expressions for the EZ depth which have the form a oc 
Ri,' n , where Ri, = 6 Z;A0/w, 2 and 1/4 < n < 1. No hysteresis occurs 
in such expressions, if we assume (see eqn . 2a) that a t , the 
inversion to surface heat flux ratio, is constant because then 
We/w, oc l/Ri,. Unlike such expressions, eqn. (22) shows 
qualitatively appropriate hysteresis. In fact, by combining 
eqns. (16), (21), (22), and (26), we see that 

a u« 4 w. 



We/w, w'9' 0 J3f 2 Zj 3 



Far from being constant, aw, /We will typically increase initially 
as surface heating and u, grow in the young day-time boundary 
layer, but will diminish later as z t becomes large. 

If we require more accuracy, we can couple eqns. (16, 18, 19, 20, 

and 26) with 



j 2uAv We av 2 a0 ( j 

H t = We + (R cb /B) | | (28) 

j f ZfA0 a 0“ j 

L J 

to provide a complete entrainment model. Here the second term in 
eqn. (28) gives the growth rate of Ah and thus a prognostic 
estimate for a. 

More physics may be added by revising the We expression to 
include the small TKE storage rate term, C st w. 2 /Ah, from eqn. (16). 
However, our present focus is amending w, , rather than refining 
the particulars of bulk entrainment models; so we choose to rest 
with the above diagnostic and prognostic formulations for a. 

In either case, hysteresis in a versus We/w, is evident, and like 
Kamada (1988b, figs. 1 and 2), depends on changes in a0 and z x . 
Both models include shear as well as heating due to entrainment, 
and avoids ad hoc surface temperature distributions, or fixed 
values of a,. Also note that eqns. (16) and (21) lets us specify 
a ( which, together with eqn. (22) , determines fl from eqn. (4) . 



BAROCLINIC SHEAR 

If the geostrophic wind varies linearly with height according to 
the thermal wind relations, u„ z = -(g/fT 0 )T y , and v g z = (g/fT 0 )T x , 
then from eqn. (23) the baroclinic stress profile must be 
parabolic, rather than linear as in the barotropic case. 

Wyngaard (1985) gave the following simple estimate for baroclinic 
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shear based on the thermal wind, 



u 



b,z 



b l U g,z 



b 2 V g,z 



/ 



V 



b,z 



biU 



*• S'Z 



b, v 

l g,z 



(29a, b) 



where b[ = 1/ (m 2 +1), b 2 = m/ (m 2 + 1), m = fZj/w,, and f is the 
Coriolis forcing. If so, the stress can be written as 



^ b ; z (gz pfT y + 2 fT p y ) 

-Uj'w' = u, 2 + — . (30) 

2 pfT 

So the baroclinic shear contribution becomes 



H gz,(b ; T x - b,T ) (b 2 gpz 2 T v + 6pTu, 2 ) 

J -u/w'u.^z = ; , 

0 6 PfT- (31) 

where terms in b, - 1 = -m 2 / (m 2 + 1) were considered neglectably 
small. The baroclinic integral is taken over the whole vertical 
range 0 to z- rather than -L to z { because the SL shear analysis 
only accounts for barotropic effects based on SL similarity. 
Scaling analysis of eqn. (31) shows that the squared and cubic 
terms in z x are small. If so, to first order, the baroclinic 
shear contribution from eqn. (31) reduces to, 



H 

1 -u.'w'u.zdz 



0 



(Z,/Z b ) w , 3 



(32) 



where we can define a baroclinic length scale as, 

fT w, 3 w, 3 

z„ . = — < 33 > 

“■■glM, - 1. 1 t v ) u.- u b I 

Usually, z b >> z r But in coastal areas mesoscale temperature and 
pressure gradients are often large. For a coastal z { of 300 m, a 
3 °C/100 km horizontal cross-wind temperature gradient, We near 
zero, f = 10‘ 4 sec' 1 , w, = 0.8 m s' 1 , and u, — 0.4 m s' 1 , then z b — 

330 m, and w, s may be enhanced baroclinically by - 50%. 
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SUMMARY 



In summary, we account for surface layer (SL), entrainment zone 
(EZ) , and baroclinic (BC) shear. Each source is given in terms 
of a dimensionless ratio of the relevant length scales z 0 , L, Ah, 
Z;, and z b . When summed and multiplied by w», this results in a 
turbulent velocity scale which includes both buoyancy and shear 
generation of TKE, namely, 



w 



*s ~ 



n - (1.25 - Ln R)L/z, + 0.5Ah/Zj + Z;/z b 

n 



- 1/3 

W, 



(34) 



Here n 



1 - a, (1 + a) 

f 

2 



a l 



We 



WeAe/w'0'o 
(C c - C d ) 

; / z ; 

Ba 0 - a u j" / a h 



the Zj/surface heat flux ratio, 



the ML entrainment rate, 



a 0 = w'0'o/We - r u Zj(l + a) 



a = Ah/z-, = 



2 . 3 (au 2 +av 2 ) 1 



r u ft zf 



au = We(u./fZj ) 2 , av = u, 2 /fZj 



R = (1 - 150z 0 /L) 1/3 - 1 



the EZ 0 jump, 

1/2 

1 

- , the EZ/BL depth ratio, 

2 



the EZ velocity jump, 



an SL shear generation term, 



w. 



3 



z 



b - 



u, 2 u 



/ 

b,z 



the baroclinic length scale. 
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